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We consider a superintegrable quantum potential in two-dimensional 
Euclidean space with a second and a third order integral of motion. The 
potential is written in terms of the fourth Painleve transcendent. We 
construct for this system a cubic algebra of integrals of motion. The algebra 
is realized in terms of parafermionic operators and we present Fock type 
representations which yield the corresponding energy spectra. We also 
discuss this potential from the point of view of higher order 
supersymmetric quantum mechanics and obtain ground state wave functions. 



1 Introduction 

Over the years many articles have been devoted to superintegrable systems 
with second order integrals of motion [1-12]. Integrable and superintegrable 
systems with third order integrals have also been studied, albeit to a lesser 
degree [13,14,15,16,17,18,19]. This article is the second in a series [18] 
devoted to superintegrable systems in quantum mechanics in 
two-dimensional Euclidean space E 2 . All classical and quantum potentials 
with one second and one third order integral of motion that separate in 
Cartesian coordinates in the two-dimensional Euclidean space were found 
by S. Gravel [16]. There are 21 quantum potentials and 8 classical ones. The 
systems investigated were of the form 
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H=^ + !f + g 1 (x)+g 2 (y) , (1.1) 

A=S.-3- + gi (x)-g 2 (y) , (1.2) 
B = ^ Aij k {L l 3 , piP2} + {h(x,y),pi} + {l 2 {x,y), P2 } , (1.3) 

i+j+fc=3 

where {, } is an anticommutator, L 3 = xP 2 — yP\ is the angular momentum. 
The constants A^ k and functions V, l± and l 2 are known [16]. 
The quantum case contains very interesting potentials written in term of 
higher transcendental functions. The irreducible potentials with rational 
functions were studied [18]. Polynomial algebras [18-27,29,30,31] and the 
parafermionic realizations of these algebras were found. The parafermionic 
realizations made it possible to construct Fock type representations and to 
obtain the energy spectra. We also studied these potentials from the point 
of view of the supersymmetric quantum mechanics [32-41]. 
Among the 21 types of superintegrable quantum potentials 5 of the 
irreducible ones are expressed in terms of Painleve transcendents [42]. Let 
us present one of the superintegrable potentials of Ref.16 written in terms 
of the fourth Painleve transcendent P±{z, a, 0) : 

uj 2 2 fku j fuj . uh 2 [uj . r— fuJ Huj 

(1.4) 

g 2 (x) = ^-y 2 , (1.5) 



where e = ±1, /' = z = J%x 

f\z) = Q^ + ^f(z) + 4zf(z) + 2(z 2 -a)f(z) + ^- ) , (1.6) 



f{z) = P 4 {z,a,f3). (1.7) 

The six Painleve transcendent functions appear in the theory of nonlinear 
differential equations. The occurence of Painleve transcendents as 
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superintegrable potentials seems somewhat surprising. It is less so once we 
remember the relation between the Schrodinger equation and the 
Korteweg-de Vries equation [43]. Solutions of the KdV include Painleve 
transcendents. Unidimensional potentials expressed in terms of Painleve 
transcendents were also obtained in the context of the dressing chains 
method [44,45,46] and conditionals and higher symmetries [47]. An 
important aspect of the fourth Painleve transcendent is the existence of 
particular solutions in terms of rational functions and classical special 
functions for very specific values of the two parameters a and (3 [48]. 
All Hamiltonians of Ref.16 are, by construction, the sum of two 
unidimensional Hamiltonians (H = H x + H y ). All the quantum potentials 
with rational function were related to supersymmetric quantum mechanics 
[19]. Higher order SUSYQM and shape invariance have been investigated 
[49,50,51,52,53,54,55]. In the case of the potential given by Eq.(1.4) and 
(1.5) the Hamiltonian H y is the well known harmonic oscillator. The 
Hamiltonian H x the corresponding Schrodinger equation has been obtained 
as a special case of third order shape invariance and solved [51]. 
This article is organized in the following way. In Section 2 we construct the 
Fock type representations for the superintegrable potential given by the 
Eq(l.l) by the means of realizations of cubic algebras in terms of a 
parafermionic algebra. In the Section 3 we will recall some aspects of third 
order shape invariance that are related to the potential given by Eq.(1.4) 
and (1.5) with e = — 1. We will also treat the case with e = 1. We will relate 
these results to those obtained using the approach involving the cubic 
algebra. In Section 4 we will consider special cases and apply results of 
Section 2 and Section 3. 



2 Cubic and parafermionic algebras 

We consider a quantum superintegrable Hamiltonian in E2 involving the 
fourth Painleve trascendent. We have two cases e=l and e=-l (with u> > 0) 



with gi(x) given in (1.4). This Hamiltonian has two integrals of motion. 
The one of the second order is given by Eq.(1.2) and Eq.(1.4). The third 
order one is given by the following equation : 

l l u 2 l h 2 u 2 

B = ~{L, Pl}+-{—x 2 y-3ygi{x), P x } -{—g lxxx {x) + {—x 2 -3gi{x))gi x {x), P y }, 

(2.2) 



where L = xP y — yPx. 

The operators A and B generate the following cubic algebra 

[A, B] = C [A, C] = 16uj 2 h 2 B (2.3) 
[B, C] = -2h 2 A 3 - 6h 2 HA 2 + 8h 2 H 3 

u 2 h* 



-(4a 2 - 20 - 6/3 - 8ea)A - 8uj 2 !fH 



8a 3 - 24a - 36a/? + 24ea 2 + 8e + 36e/5) 



The Casimir operator can be written as a polynomial in the Hamiltonian 

K = -16h 2 H 4 + ^——{4a 2 -8a + 4-a(3)H 2 (2.4) 
3 
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8a 3 - 24ea 2 + 24a + 36a/5 - 8e - 36e(3)H 
(4a - 8ea - 8 - 6/3) . 



Realizations of cubic algebras in terms of parafermionic algebras have been 
discussed in our previous article [19]. Our potential belong to the Case 2 of 
Ref.19. The cubic algebra has the form : 

[A, B] = C, [A, C] = SB, [B, C] = fiA 3 + vA 2 + £A + ( , (2.5) 



where 

H = Ho, v=v + viH, C = + ZiH + &H 2 (2.6) 
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C = Co + CiH + C2H 2 + C3H 3 , 6 = 5 + 5 1 H 



This algebra has been realized in terms of a deformed oscillator algebra of 
the form 

[N,b t ] = b\ [N,b] = -b, b t b = <S>(N), 66' = $(iV + l) . (2.7) 
The structure function $(iV) is given by 

H=¥ + { + + nf + C-f-f)(N + „)« + (f)(N + u)* . 

We can use Eq.(2.4) to rewrite the structure function in terms of the 
Hamiltonian. 

2.1 Case e = 1. 

From the general formula we obtain for our particular case the following 
structure function for e = 1 

•(*) = -4^(* + u - ( + + „ - + I - |)) (2.9) 

+ " - + 5<° + 2 - 3i /f »><* + " - < 2^ + S< a + 2 + 3! /f »> ' 

To obtain unitary representations [25,26] we should impose three 
constraints given by 

$(p + l,Ui,fc) = 0, $(0,u,Jfe) = 0, 0(x) > 0, V x>0 . (2.10) 

We have to distinguish the two cases j3 < and /? > 0. For (3 < we get 
four possible values for u with $(0, u,k) = 



-E 5 a -El, -0. 
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-E \, M3, E 1 

M3= 2^ + 6 (a + 2 - 3 V^ ) ' " 4= 2^ + 2 

We insert all these solutions for u and apply the constraint 
$(p + 1, Ui, k) = 0, with i=l,2,3,4 to find the energy spectrum. 
Case 1 



= 4hWx(p + 1 - x){x+ 1 --^-^ ){x+ 1 --^ + ^ ) (2.12) 

4 a 

E = hLu(p+---). (2.13) 

Case 2 



^MP+g+ 6 + V 8 } - (2 - 15) 



Case 3 



^=^(p+| + f-V-r ) - (2 - 17) 



Case 4 

We get three solutions for this case with negative energy 



1 N , a , /-/? , 1 N , a /-/? 



$(x) = 4liVz(p+l-x)((p+-)--- t /^-x)((p+l)-^-x) ) (2.20) 



$(a;) =4^ 4 o; 2 x(p + l-x)((p+-) + - + y — -) + - - ^ — -x), 

(2.18) 

£=-fiu;(P + ^ + f), (2.19) 
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E = -hu(p+ 7 --^- ] ff), (2.21) 



<S>( x )=4hWx(p+l-x)((p+^)~ + ^-x)((p+l) + ^-x), (2.22) 

^"^V^/t 1 ' (2 - 23) 



To obtain unitary representions we should also impose <j){x) to be a real 
function and (f)(x) > for x > 0. The constraints do not allow all values for 
a and /? so it may happen that only some of the states are physically 
meaningful. We can have 1, 2 or 3 infinite sequences of energies that 
correspond to each unitary representation. 
For f3 > we have two solutions 

8 n 2 <y 1 

$(x) = AhWx{p + 1 - x){x 2 + (1 - a)x - - + ^ - - + =■), (2.24) 

8 4 2 4 

£ = MP+^-f), (2-25) 

8 a 2 a 1 

$(x) = 4/i 4 a; 2 2;(p + l-2;)(x 2 -(l + a + 2p)x + p 2 + ap + p-^ + — + - + -), 

(2.26) 

2 a 

£=-fcj(p + - + -). (2.27) 



2.2 Casee=-1. 

For the case e=-l we obtain the following expression for the structure 
function 

E 1 — E 1 a 

•(*) = + « - (— + -))(x +„-(+)) (2.28) 
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E=-hu(p+- + -), 



(2.37) 



$(x) = 4^ 4 o; 2 a;(p+l-x)((p+-)----\/^ — x)((p+l)-\/ — x), (2.38) 



tf=-MP+j|-£-V-/), (2-39) 



$(x) = 4^ 4 o; 2 a;(p+l-x)((p+-)---'y/^--x)((p+l) + -y/^f -x), (2.40) 



^"^V^ + VT 1 ' (2 - 41) 

One interesting aspect of this potential is that we can have three, two or 
one series of equidistant energy levels. 
For (3 > we get the following solution 

$(x) = 4ftVx(p + 1 - x)(x 2 - (1 + a)x - - + — + - + h, (2.42) 

8 4 2 4 

2 a 

E = hw(p+---), (2-43) 

/? a 2 a 9 

$(x) = 4^, 4 a; 2 x(p + 1-x) (x 2 - (3 + a + 2p)x + p 2 + ap + 3p- ^ + X + "2 + 4^' 

(2.44) 

£ = _^(p + l + |). (2.45) 



3 Third order shape invariance and 
superintegrable systems 

The concept of higher-derivative supersymmetric quantum mechanics 
(HSQM) was introduced by A.A.Andrianov, M.V.Ioffe and V.P.Spiridonov 
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[49]. HSQM is characterized by polynomial relations between supercharges 
and the Hamiltonian. Second order derivative supersymmetry was 
investigated in the Ref.50. We will present in this section results not given 
in Ref.51 but directly related to the potential given by Eq.(l.l) with e — 1. 
Let us recall some aspects of the particular case of third order shape 
invariance related to the potential with e = — 1 obtained in the Ref.51. In 
SUSYQM two superpartners are isospectral or almost isospectral and if we 
know the spectrum and the eigenfunctions of one superpartner we can 
obtain the spectrum and the eigenfunctions of the other superpartner. A 
special case occurs when the two superpartners Vi(x, ao) and V2(x,ao) 
satisfy the relation V 2 (x, ai) = Vi(x, a ) + -R(ai) where ai = /(ao) and R(ai) 
do not depend on x. In this special case we can find directly the energy and 
the eigenfunctions. The superpartners are called shape invariant potentials 
(SIP). We consider the following particular case of shape invariance 

H x a} = a! (Hi + 2A) , (3.1) 

where a) and a are third order operators. This particular case of shape 
invariance can be constructed from a first order and second order 
supersymmetry given by the following interwining relations 

HiJ = q\H 2 + 2X), E 1 M^ = M^E 2 , (3.2) 

where 

H t = P* + Vi(x) , (3.3) 

qi = d + W(x), q = -d + W(x) , (3.4) 

M f = d 2 -2h(x)d + b(x), M = d 2 + 2h(x)d + b(x) . (3.5) 

The key element in obtaining the equivalence between Eq.(3.1) and Eq.(3.2) 
is to define the following third order operators a and a) written as products 
of first order and second order supercharges 

a) = q j M, a = M ] q . (3.6) 

The third order shape invariance of the form given by Eq.(3.1) can be 
investigated using Eq.(3.2). The two interwining relations of Eq.(3.2) give 
respectively the following relations 

Vi = W'(x) + W 2 (x), V 2 = -W'(x) + W 2 (x) - 2A (3.7) 
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and 

^ = T2ft ' W + ^) + ^£l-|M__l_ +7 , (3 . 8 ) 

+ + ^ + . (3 . 9) 

Eq.(3.7), (3.8) and (3.9) impose that the potential V\ should have the form 
Vx = -2h'(x) + 4h 2 (x) + AXxh{x) + X 2 x 2 - X , (3.10) 

with 

h"(x) = f ^ + 6h 3 (x) + 8Xxh 2 (x) + 2(X 2 x 2 -(X + 1 ))h(x) + -^, (3.11) 
2n{x) 2h[x) 

W(x) = W 3 (x) = -2h(x) - Xx . (3.12) 
As in the case of first order supersymmetry we can define 

«-(£?) f-Qfj • w 

We get the following SUSY-algebra 

[H,Q] = [H,Q ] ]=0, {Q,Q} = {Q\Q ] }=0, {Q,Qi} = (H- 1 ) 2 + d . 

(3.14) 

Eq.(3.11) can be transformed into the equation for the fourth Painleve 
transcendent (1.4) by the following transformations 

h(x) = ±V\f(z), z = V\~x, « = 1 + I P=™, A = |, (3.15) 
and we obtain 

(3.16) 

Vi(x) is the x part of the potential in (1.1) and coincides with gi(x) in 
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Eq.(1.4) up to a constant. A particular case of third order shape invariance 
called « reducible »was considered in Ref.51 by imposing further conditions. 
These conditions are d < and the existence of real functions W\ and W 2 
such that 

= (d + W 1 (x))(d + W 2 (x)), W lt2 = -h{x) ± h ^ ~ , (3.17) 

(reducible means that M * factorizes into product of two first order 
operators with real functions). The spectrum was obtained for cases where 
normalizable zero modes of the annihilation operator exist. Zero modes of 
the annihilation operator satisfy 

a^ 0) =0 , (3.18) 

(we use the terminology of HSQM where « zero mode » refers to Eq.(3.18) 
so that zero modes may not have energy Eq — 0). The energies of zero 
modes were obtained by imposing the vanishing of the norm of aip^ which 
involves the average of the operator product a^a. 

a}a = q^MM ] q = q\(H 2 - 1 f + d)q = H 1 ((H 1 - 1 -2Xf + d) . (3.19) 
The energies of the zero modes are 

e{ 0) = 0, Ef ] = 7 + 2A + V^, E^ = 7 + 2A - V^d . (3.20) 
The corresponding eigenfunctions ipf^ can be calculated explicitly and are 

^0( x )= e rW<fe' ; (3 _2i) 

$(x) = (W 2 (x) - W 3 (x))e~f XW2{x ' )dx ' , (3.22) 
i/j° 3 (x) = (2V^+(W 2 (x)-W 3 (x))(W 1 (x) + W 2 (x)))e-^ m{x ' )dx ' . (3.23) 



The creation operator can also have zero modes (p k which correspond to a 
possible truncation of the sequence of excited levels. They were obtained by 
considering the following product 

aa) = (Hi +2A)((#i - -jf + d) . (3.24) 
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The energies of the zero modes are 



E$»=-y-V=d t E$ s >=-y + V=d t E^ = -2X 



(3.25) 



with the corresponding eigenfunctions 



$(x) = (W^x) + W 2 {x))e^ W2{x ' )d: 



(3.26) 



(3.27) 



= h + 2\ + ^+(W 1 (x) + W 2 (x))(W 2 (x)-W 3 (x)))e-f Xw ^ dx ' . 

(3.28) 



For non singular potentials it is not possible to have the negative energy 
and the total number of zero modes of the annihilation and creation 
operator cannot be more than three because of the asymptotics of the 
eigenfunctions. We can have three, two or one infinite sequence of levels. 
These results coincide with those obtained as from the analysis of Fock type 
representations of the cubic algebra of the superintegrable potential. When 
we apply the creation operator on zero modes we create eigenfunctions 
with 2A more energy. These energies are corroborated (when we add a 
harmonic oscillator in the y direction) by those obtained using the cubic 
algebra and given by Eq.(2.31), Eq.(2.33) and Eq.(2.35). 
When a potential allowes only one infinite sequence of energies, this 
potential may also allow a singlet state or a doublet of states 



From an algebraic point of view these states correspond to trivial 
irreducible representations. Such a case was discussed in Ref. [19] for the 
potential V = h 2 ( x ' ^ + ^ x \y + ( x + a y ) ■ This potential is a special case of 
the potential given by the Eq(l.l). The observed singlet state can now be 
naturally understood as a phenomenon of third order shape invariance. 
All the results we presented apply to our potential for e = —1. We will 
present here the results that will be applicable to the case e — 1. We follow 
the same approach as in Ref. 51 and we consider the following potential 



The Eq.(3.9), Eq.(3.11) and Eq.(3.12) remain the same. We can define as 
for Vi in Eq.(3.16) a potential V 2 using the transformations of Eq.(3.15). 



a + ip(x) = a ip(x) = 0, (a + ) 2 ip(x) = a ip(x) = 



(3.29) 



V 2 = 2h'{x) + Ah 2 {x) + A\xh(x) + X 2 x 2 - A 



(3.30) 
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The Hamiltonian H 2 of the form given by Eq.(3.3) with the potential given 
by Eq.(3.30) satisfies 

H 2 a) = a\H 2 + 2\) . (3.31) 

when we postulate 

a) = Mq\ a = qM* . (3.32) 
We have the following product 

a)a = H 2 ((H 2 -j) 2 + d), aa f = (H 2 + 2X)((H 2 + 2A - 7 ) 2 + d) .(3.33) 

The energies of the zero modes of the creation and annihilation operator 
are obtained by imposing the vanishing of their norm. This involves the 
average of the operator product a) a and ao) given by Eq.(3.33). The 
eigenfunctions of the zero modes are obtained directly by solving aipf^ = 
and = 0. The energies of zero modes of the annihilation operator are 

£j 0) = 0, ^= 7 -V=d, 4^=-r + V=d , (3.34) 
with the corresponding eigenfunctions 

$>(x) = (^-V^+(W 1 (x) + W 2 (x))(W 1 (x)-W 3 (x)))e^ w ^ dx ' , (3.35) 

$(x) = (W^x) + W 2 (x))e-^ Wl{x ' )dx ' , (3.36) 
1%( x ) = e-S mw *W . (3.37) 
The energies of zero modes of the creation operator are 

E{ 0) = -2A, E [ 2 Q) = 7 - 2A - V^d, Ei 0) = 7 - 2A + , (3.38) 

with the corresponding eigenfunctions 

(t P i{x) = e -Fw z{x > )dx > ; (3 39) 

$(x) = {W l {x) - W 3 (x))e^ Wlix ' )dx ' , (3.40) 

$j(x) = (-2V^+(W 1 (x)-W 3 (x))(W 1 (x) + W 2 (x)))e^ W2 ^ dx> . (3.41) 

Again the total number of zero modes of the annihilation and creation 
operator cannot be more than three because of the asymptotics of the 
eigenfunctions. We can have three, two or one infinite sequence of levels. 
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When we apply the creation operator a) on zero modes we create 
eigenfunctions with 2A more energy. These energies are corroborated (when 
we add a harmonic oscillator in the y direction) by those obtained by the 
cubic algebra and given by Eq.(2.13), Eq.(2.15) and Eq.(2.17). When a 
potential possess only one infinite sequence of energies, this potential may 
also possess a singlet state or a doublet states. 

We will discuss the irreducible case that appears when d > 0. For Vi(x) we 
get 

E®=0, ^° (x) = e^ w ^ dx '. (3.42) 

For V2(x) we get 

E[ 0) =0, ip°{x) = (^-\Q+(W 1 (x) + W 2 (x))(W 1 (x)-W 3 (x)))e^ w ^ x ' )dx '. 

(3.43) 

and W\ and W 2 are not real functions. 



4 Special cases 

The fourth Painleve transcendent satisfying Eq.(1.6) depends on two 
parameters and special solutions in terms of rational or classical special 
functions exist [48]. In this section, we will give the unitary representations, 
the corresponding energy spectra and the eigenfunctions for some special 
cases. 



4.1 Case a = 5, f3 = -8, f(z) = A \^ + ^ z l^ and e = 1. 
We have with Eq.(1.4) and (1.5) 



2 



8ti 3 u Ah 2 u 2hu 



From the cubic algebra we get two unitary representations. The first 
unitary representation is given by Eq.(2.14) with the corresponding energy 
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given by Eq.(2.15) 

g 

^(i) = 4fi 4 w 2 i(p+l-i)(i + 3)(i + 2) 1 E = hu(p+-), (4.2) 

The second solution is given by Eq(2.12) or Eq.(2.20) with the 
corresponding energy spectrum 

<j)(x) = 4hWx(p + 1 - x) (x - 3) (x - 1) , E = fajj{p — -). (4.3) 
This representation is valid only for p=0. 

We will also treat this systems using the results on supersymmetry. The 
eigenfunctions for the x part consist of an infinite sequence xj) n (x) starting 
from psi\{x) of Eq.(3.37) and a singlet state x( x ) given by Eq.(3.35) and 
(3.40) 

ax(x) = 0, a^x(x) = . (4.5) 

The creation and annihilation operator are given by Eq.(3.32) with the 
following expressions for W±, W 2 and W 3 

-(-h + 2ux 2 )(9h 3 + 27h 2 cox 2 + 12Huj 2 x 4 + 4oAr 6 ) 
1_ hx{3h + 2ux 2 ){3h 2 + 4u 2 x 4 ) ' 

-ifi - 2ux 2 )(3h 2 + 3hux 2 + 2u 2 x 4 ) . s 

Wo = — — (4 7) 

2 hx(3h 2 + 8hiux 2 + Au 2 x A ) ' y ' 1 

-ux(-9h 3 + 22h 2 cux 2 + 20frw 2 x 4 + 8u 3 x 6 ) 
3 ~ h(h + 2ujx 2 )(3h 2 + Au 2 x A ) ' { ' 

With the eigenfunctions for the y part of the Hamiltonian and the formula 
for the energies given by Eq.(3.34) we obtain the two following series of 
solutions 

ifJn,k = ^n{x)e~^H k (J l ^y), E = hw(n + k + -) , (4.9) 
K = x{x)e~^H m (J^y), E m = ?iw(m--). (4.10) 



16 



4.2 Case a = 5, j3 = -8, /(z) = ^p+ipp+j^ and e = -1. 

We have with Eq.(1.4) and (1.5) 

, , c^ 2 2 , 192fr 4 ^V 16fr 2 u; 2 x 2 

V) = + V ) " (i^i + 3fi2)2 + 4^ + 3^ ' (4 " n) 

From the cubic algebra we obtain three unitary representations. The first 
unitary representation is given by Eq.(2.32) with the corresponding energy 
spectrum Eq.(2.33) 

4>{x) =4hWx{p+ 1 -x)0 + 4)(x + 2), E = frw{p + 3), (4.12) 



The second solution is given by Eq(2.38) with the corresponding energy 
spectrum given by Eq.(2.39) 

(j){x) = Ah A u 2 x(p + 1 - x)(p - 3 - x)(p - 1 - x), E = -huj(p - 1). (4.13) 
This representation is valid only for p=0,l. 

The third solution is given by Eq(2.30) with the corresponding energy 
spectrum given by Eq.(2.31) 

^(i) = 4ftVi(p+l-i)(i-3)(i-2), E = kuj(p-l). (4.14) 
This representation is valid only for p=0,l. 

We investigate this system using the results on supersymmetry. The 
eigenfunctions for the x part consist of an infinite sequence tjj n (x) starting 
from ^o(x) given by Eq.(3.22) and doublet states Xi( x ) an d X%{ x ) given by 
Eq.(3.22), (3.21) and (3.26) 



axi(x) = 0, a'xi(x) = Xnfa), a^fa) = . (4.17) 
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The creation and annihilation operators are given by Eq.(3.6) with Wi, W 2 
and W 3 as in Eq.(4.5), Eq.(4.6) and Eq.(4.7). 

With the eigenfunctions in the y part and the formula for the energies given 
by Eq.(3.20) we obtain the three following kinds of solutions 

i>n,k = Mx)e~^ H k(J^y), E = hw(n + k + 3) , (4.18) 

<f>mi = Xi{x)e~^ H mi (J —y), E mi = fuj(mi - 1) , (4.19) 

4>m 2 = X2(x)e~^ H m2 (^^y), E m2 = hum 2 . (4.20) 

4.3 Case a = 0, (3 = -§, f(z) = —\z and e = 1. 

We have 

V(x,y) = ^x 2 + y 2 ) (4.21) 

From the cubic algebra we get three cases with unitary representations and 
using Eq.(2.12) to Eq.(2.17) we have 

(f)(x) =4h 4 u> 2 x(p + l -x)(x + h(x+ ~), E = hw(p+^-), (4.22) 

<p(x) =4h 4 u 2 x(p + l-x)(x--)(x + -), E = huj(p+1), (4.23) 

3 3 

2 1 2 

<j>{x) =4hWx(p+l-x)(x- ~)(x- -), E = huj(p+ -). (4.24) 

3 3 3 

We apply the results coming from supersymmetry. we get the following 
known eigenfunctions from Eq.(3.35), Eq.(3.37) and E.(3.36) respectively 
and the corresponding energy with the Eq.(3.34) 



VWi = N nikl (aT 1 e—(-3h + 2ux 2 )e~H k U-y), (4.25) 

4 

E 1 = hu(n 1 + k 1 + -), (4.26) 
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-ujx 2 uiy 2 jUJ 



</W 2 = N n2k2 (a^e^xe-^rH k2 (^y), (4.27) 

E 2 = hu(n 2 + k 2 + l). (4.28) 

ipn 3 ,k 3 = N n3k3 (a^e^e-^H k3 (^y), (4.29) 

E 2 = hu(n 2 + k 2 + -), (4.30) 

The creation and annihilation operator are given by Eq.(3.32) with the 
following expression for W\, W 2 and W3 

1 LUX TTr 1 UJX TTr LUX 

Wi = ~ + —, W 2 = — + —, W 3 = -—. 4.31 
x in x in in 

4A a = -l,/? = -f, f(z) = -f-^fe|) and c = 1. 
We have 

u 2 A „ 9 , 24ft 3 u; 4fr 2 w 

2 9 [2ojx z + 3rt) 2 (2u;x^ + in) 



From the cubic algebra we get the three cases with unitary representations 
from Eq.(2.12) to Eq.(2.17) 

0(x)=4/iVx(p + l-x)(x+-)(x+-), £ = Mp+~), (4.33) 

14 4 
<f>{x) = 4ftVx(p + 1 -x)(x - -)(x + -), £? = /jw(p+-), (4.34) 

5 4 

0(x) = 4h 4 u 2 x(p + 1 - x)(x - -)(x - -), E = huo(p). (4.35) 



From the supersymmery we obtain the following eigenfunctions from 
Eq.(3.35), Eq.(3.36) and E.(3.37) respectively and the energy with the 
Eq.(3.34) 

, +vn 2 (-45fr 2 + 4u; 2 x 4 ) W y 2 tt . /o7 . 
VVu,* = iVnA^)"^— x l (3n+2u;x2) J e- — H kl (^-y), (4.36) 
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E 1 = huj(n 1 + k 1 + ~), (4.37) 



. Uno (9h 2 - 12hux 2 - 4a; V) ^ TT . fu . fl . 

*.* = A^(.»r«-ri (4.38) 

4 

£ 2 = ^j(n 2 + A; 2 + -). (4.39) 



^ = ^(a^-p^ye^F^^,). (4.40) 

E 3 = hu(n 3 + k 3 ), (4.41) 

The creation and annihilation operators are given by Eq.(3.32) with the 
following expressions for Wi, W2 and W 3 

-27K 3 + 27h 2 ux 2 + A8huj 2 x A + Au 3 x 6 
1 ~ 27ft 3 x - 3QH 2 ujx 3 - 12Hu 2 x 5 ' ^ ' ' 

_ 351h 3 ux + 126h 2 u 2 x 3 + 12huj 3 x 5 - Stu 4 x 7 
2 ~ 81H 4 - 54h 3 ux 2 - W6h 2 LU 2 x 4 - 24hu 3 x 6 ' 
-9h 2 - 3hux 2 + 2u 2 x 4 
Ws = 9h 2 x + 6hux 3 ' ^ AM > 

4.5 Case a = 0, /? = -2, /(z) = -2z - ^(z) and e = 1. 
We have 

*W = ^, #) = 1-^) • (4.45) 
E c (z) is the complementary error function and is given by 



00 

2 



E c (z) = 4= / e _ta dt (4.46) 



/7T 

We have with Eq.(1.4) and (1.5) 

a; , 9 , N 2,_ 4e * tnw 4e >> + u;V'^£ 

y) = — (x 2 + y 2 ) - -fej + — — . = NN9 + 



2 3 7r(l-^ c (yfx)) 2 7r(l-tE c (yfx)) 

(4.47) 



20 



The cubic algebra provides two unitary representations. The first unitary 
representation is given by the Eq.(2.12) with the corresponding energy 
given by Eq.(2.13). The Eq.(2.14) and (2.15) give the same unitary 
representation and energy spectrum and we have 

4 

<f>(x) = Ah A uj 2 x(p+ 1 - x)x(x + 1), E = hu(p+-), (4.48) 

3 

The second solution is given by Eq(2.16) with the corresponding energy 
spectrum given by Eq.(2.17) 

<j)(x) = AhWx(p+ 1 -x)(x - l) 2 , E=tuj(p+-). (4.49) 
This representation is valid for p=0. 

We also use supersymmetry to treat this system. The eigenfunctions for the 
x part consist of an infinite sequence ifj n (x) starting from if)® given by 
Eq.(3.37) and a singlet state x{ x ) given by Eq.(3.36) and (3.39) 



, , x Ar / +\n -3^ 2 (-Wfiio - e n v^Fwx + e » s/ntuxE c (^ffx)) 



u(-l + tE c (y/%x)) 



(4.50) 



= - ' <451) 

a X {x) = 0, a)x{x) = . (4.52) 

The creation and annihilation operators are given by Eq.(3.32) with the 
following expressions for W±, W2 and W3 



(-tux - e^^/7i^(h + lux 2 ) + e^- v / 7rt A /f(^ + u)x 2 )E c (y/%x)) 
h(-t-e^f^^/fx + e^ v ^+ ^/%xE c (^%x)) 

(4.53) 

W 2 = ^ , (4.54) 
A = (2t 2 VujE + 3e"^~ y/ntux - e ^ 7l \j^( h ~ ^ 2 )+ (4-55) 
(-3e~ y/jtfux + 2e~7rtJ -(h - ux 2 ))E c (J -x) 
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1—1 1 LUX 2 ^ FtO U,X 2 _ /U7 . /U7 s . 

B = h v ^(-l+tE c (J-x)){-t-e s ^J-x + e * y/ntJ -xE c (J -x)), 

(4.56) 



2e fi + 

W 3 = — + v ^ . (4.57) 

/l ^ 0F(1 - t£ c (,/f x ) 



With the eigenfunctions for the y part of the Hamiltonian and the formula 
for the energies given by Eq.(3.34) we obtain the two following families of 
solutions 

^n,k = ^n(x)e-^H k (^y), E = hw(n + k + ^) , (4.58) 
m = x(i)e"^ff m (J^y), £ m = foj(m + -). (4.59) 



4.6 Case a = 0, /? = -2, /(z) = -2z - #(z) and e = -1. 

This case give the harmonic oscillator. The zero mode is given by Eq.(3.22) 
is the well known ground state of the harmonic oscillator. There is other 
special solutions in terms of the complementary error function exist. 
Many special solutions of the fourth Painleve equation will give us singular 
Hamiltonians. These potentials can be regularized in several manners 
[19,56,57,58]. 



5 Conclusion 

The main results of this article are that we have constructed the cubic 
algebra, Fock type presentations and the corresponding energy spectrum for 
the potential given by Eq.(1.4), (1.5). Other superintegrable potentials 
written in term of Painleve transcendents are known [16] namely : 

V 1 = tffalPjfax) + uJ 2 2 Pi(u 2 y)) (5.1) 
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V 2 = ay + h 2 u?Pi{u)x) (5.2) 
V 3 = bx + ay + (2hb)lpfj((^x, 0) (5.3) 

V A = ay + (2hW)HPU(^x, k) + P£((^)H k). (5.4) 

These potentials together with that in Eq.(1.4) and Eq.(1.5) were also 
obtained as one dimensional potentials in the context of higher and 
conditional symmetries by W.I.Fushchych and A.G.Nikitin [46]. For these 
four superintegrable potentials the simplest underlying structure of the type 
(2.5) is actually a finite dimensional Lie algebra that does not allow us to 
find the energy spectrum. Let us present these algebras : 
For Vi we have : 

[A, B]=C = -iffuful [A, C] = 0, [B, C] = . (5.5) 



For V2 we have : 

[A, B] = C = -ih 5 u 5 , [A,C]=0, [B,C] = . (5.6) 

For V3 and V4 we have : 

[A,B] = C = 4abih(H+-A), [A,C}=0, [B, C] = 8a 2 b 2 h 2 (H + -A) . 

(5.7) 

These algebras coincide with the classical Poisson algebras presented earlier 
[17]. In these four cases we have a triplet of commuting operators. The x 
part of the potential V4 was also obtained in the context of supersymmetric 
quantum mechanics [55]. The methods of this article are not directly 
applicable in that case, but it may be possible to generalize them. 
The question of how these aspects of SUSYQM, shape invariance and 
superintegrability are related is interesting and will require more study. 
Supersymmetry could also be a tool for the classification of superintegrable 
potentials. Higher order supersymmetry could be a suitable approach to 
treat these potentials. The search for superintegrable systems with higher 
order integrals of motion is thus closely related to the subject of polynomial 
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algebras and higher order supersymmetric quantum mechanics. 
The search for a grand unifying theory in particle physics is an important 
problem of comtemporary physics. One model that is envisaged as a 
candidate is string theory. The x part of the potential given by Eq.(l.l) 
appears also in the context of string theory [59] where supersymmetry is 
used as a method for constructing exact solutions. A more recent article 
[60] discusses how supersymmetric quantum mechanics can be used to 
construct solutions in string theory. 
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